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3-d lattice SU(3) free energy to four loops * 
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We report on the perturbative computation of the 3d lattice Yang-Mills free energy to four loops by means of 
Numerical Stochastic Perturbation Theory. The known first and second orders have been correctly reproduced; 
the third and fourth order coefficients are new results and the known logarithmic IR divergence in the fourth 
order has been correctly identified. Progress is being made in switching to the gluon mass IR regularization and 
the related inclusion of the Faddeev-Popov determinant. 



1. Introduction 

The main goal of finite temperatme QCD is 
to characterize the deconfinement transition bet- 
ween the low temperature regime, ruled by con- 
finement, and the high temperature regime ruled 
by asymptotic freedom. In this respect the free 
energy density, i.e. the pressure, is a sort of 
'theoretical laboratory' in which to study such a 
transition, and, eventually, a potential candidate 
observable in heavy ion collisions. Evidently, the 
large T limit of the pressure is that of an ideal 
gas of non-interacting particles, p cx T^; unfortu- 
nately, lattice simulations across the whole tem- 
perature range are not feasible and our efforts 
and understanding are confined to the opposite 
ends of the temperature range: namely, the low 
temperature regime (up to some 4 5 times the 
transition temperature Tc ^ 200 MeV) which is 
accessible to actual computational resources, and 
the high temperature regime at which the pres- 
sure exhibits a purely perturbative behaviour. 

Dimensional reduction has been applied [2] 
to fill the gap: Ad QCD is matched to M SU{3) 
Yang-Mills (YM) coupled to a Higgs field in the 
adjoint representation; such a theory can then be 
matched to 3d SU{3) YM, which captures the 
ultrasoft degrees of freedom; both these reduc- 
tions have been successfully performed in the MS 



scheme. Finally, 3d YM has to be treated non- 
perturbatively, the only effective means being the 
lattice. Indeed lattice computations can be con- 
sistently incorporated into the dimensional reduc- 
tion strategy: the theory is superrenormalizable 
and all divergences can be computed perturba- 
tively; this allows a clean matching of the schemes 
in the continuum. It is well known that comput- 
ing at high orders in LPT is quite hard a task; in 
this respect the approach of Numerical Stochastic 
Perturbation Theory (NSPT) ,3^ is a very efficient 
one. 

Let f = —V\nZ he the free energy associated 
with the Wilson action Sw = PoJ2pi^ ~ Hp), 
Po = 2Nc/ {a'^~'^gQ) being the dimensionless bare 
lattice inverse coupling in d dimensions. To com- 
pute the free energy one can revert to the com- 
putation of the plaquette, being 
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then, a weak-coupling expansion of the plaquette 
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uniquely determines the corresponding expansion 
for /, and for A^c = 3 and d = 3 
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so that the first four coefficients uniquely deter- 
mine the matching to the continuum; ci and C2 
being long known, C3 and C4 were the aims of 
the computation. The presence has already been 
established |1] of a logarithmic IR divergence at 
fourth loop order; it was our aim to recover the 
scheme independent coefficient of such a loga- 
rithm: the tilde on C4 reminds us that the IR 
divergence has to be isolated and subtracted by 
means of a definite IR regulator^ which we choose 
to be the finite lattice size L inherent to any lat- 
tice simulation. 

2. Computational setup 

We have inserted the expansion U^{n) = 1 + 

Si=i ^ ^M*^ ("■) Langevin equation 

dtU^ ^ [-iVS[U^] - i^p^ , (4) 

which implements the strategy of Stochastic 
Quantization on lattice gauge fields (77 is a gaus- 
sian noise). We have integrated such an equa- 
tion in the Euler scheme and linearly extrapo- 
lated to vanishing stochastic time step. We have 
computed the coefficients of the expansion (3) on 
lattice sizes ranging from L — 5 to L = 16; up to 
three loops we also performed the computations 
on a L = 18 lattice. 

We finally extrapolated to infinite lattice size 
according to the analysis by Symanzik which en- 
tails subleading logarithms, suppressed by in- 
verse powers of L; recall that the finite volume 
acts as the IR regulator needed at 4 loops and en- 
ters a leading logarithm. The indeterminacies we 
quote are dominated by this extrapolation proce- 
dure: our data do not allow a clear cut between a 
logarithm and a constant and we have been forced 
to effective extrapolations based on pure power- 
like fits; the spreading of the results (figures) cor- 
responds to different set of powers included in the 
fit. More details can be found in 

Most of the computer simulations have been 
performed on a PC cluster of 10 bi-processors 
Athlon MP2200. A programming environment 
for LGT-NSPT has been set up in C-f-l- by V. 
Miccio with massive use of (CH — h specific) classes 
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Figure 1. Extrapolating fits for C3. 



and methods well suited to handle lattice and al- 
gebraic structures. Some statistics came from a 
14 Intel Xeon 2.0 Ghz PC cluster which has be- 
come available along the way. 

3. Results 

In we present the results we obtained for the 
coefficients c^^'' at various L and the values ex- 
trapolated to infinite volume; the two benchmark 
coefficient are shown to be in excellent agreement 
with the known results. Here we simply plot (fig- 
ures 1,2) the extrapolating fits for the third and 
fourth coefficients, the new results, and quote the 
numbers: 

C3 = -Q-QO[ll-) , C4 = -25.8(4) . 

In we discuss at length about the IR di- 
vergence at four-loop order; let us just recall 
that it is quite hard to recognize the presence 
of a logarithm out of any set of data; we found 
that fits with no log had slightly but system- 
atically greater than fits with a logarithmic 
divergence, and, most of all, that we could fit 
quite well (04"^ = 1.1(2)) the coefficient 04"^ = 
81 (688 - 1577rV4)/(47r)'' = 0.9765 which has 
been computed in 0]. 
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Figure 2. Extrapolating fits for the log- 
subtracted C4. 



4. Perspectives 

A definite IR regulator enables one to isolate 
a finite part out of the 4 loops contribution to 
the plaquette; we have been employing the finite 
volume inherent to any computer simulation. A 
mass is a commonly used regulator in continuum 
computations; the coupling itself (dimensionful in 
3d YM) is quite a natural IR regulator for non- 
perturbative lattice simulations. Despite the fact 
that the coefficient of the logarithmic divergence 
is universal it is clear that definite coefficients re- 
late different regularizations. We are interested 
in the perturbative matching between the lattice 
and a continuum perturbative scheme: in both 
approaches the same gluon mass IR regulator can 
be employed so that the same mismatch shows 
up with respect to data coming from nonpertur- 
bative computer simulations with the coupling as 
IR regulator. Employing massive gluons is the na- 
tural extension of our computational programme. 

Inserting a mass into the gluon propagator 
amounts to switching to the Lie algebra for the 
definition of the new lattice action and this we 
perform at a definite step of the computation. It 
is well known that the Lie algebra YM action is 
made out of a gauge fixing term and the ghost 
action resulting from the FP determinant. Our 



NSPT strategy entails the construction of the Eu- 
ler scheme solution U' = exp{—eF)U to Eq. (4); 
the quantity F belongs to the Lie algebra and our 
algorithm entails going back and forth from the 
group to the algebra: at that stage we include 
the new contribution to the action. The strategy 
for a determinant is well established |()I7| and is 
employed by the Parma group in the unquenched 
NSPT approach to lattice QCD |S|; the determi- 
nant is simulated without introducing fermion or 
ghost fields and the inversions of the nonlocal ma- 
trix is performed perturbatively via an efficient 
FFT by a back and forth procedure from Fourier 
space: the zeroth order of any propagator is di- 
agonal in momentum space and the construction 
of successive orders does not require any matrix 
inversion. The implementation of such a strategy 
is in progress; signals of the stabilization of gauge 
dependent quantities as well as of the invariance 
of gauge invariant quantities make us confident 
about the reliability of the whole procedure. 
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